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Preface 


The  report  entitled  “Vibration  Damping  by  a  Nearly  Continuous  Distribution  of 
Nearly  Undamped  Oscillators,”  NSWCCD-70-TR-1999-120  (April  1999)  was  submitted 
as  a  paper  to  be  published  in  the  Journal  of  Sound  and  Vibration.  Two  reviewers  set  upon 
this  paper  and  I  was  asked  to  shorten  and  revise  some  of  the  arguments.  Unfortunately, 
the  reviewers  failed  to  understand  the  gist  of  the  paper,  in  part  I  suppose  the  blame  is  mine. 
Consequently,  I  decided  to  reiterate  a  major  argument  in  the  submitted  paper  and  make  it 
the  center  piece  of  the  revised  version.  The  letter  I  sent  with  the  revised  manuscript  to  the 
Chief  American  Editor  of  this  Journal;  Professor  Soedel,  is  enclosed  for  your  perusal.  I 
thought  some  readers  of  the  present  report  may  enjoy  a  reprint  of  this  revised  version  of  the 
manuscript. 


Dear  Professor  Soedel: 

Enclosed  is  a  condensed  and  revised  version  of  the  paper  previously  submitted  for 
publication  under  the  title  “Vibration  Damping  by  a  Nearly  Continuous  Distribution  of 
Nearly  Undamped  Oscillators.”  Since  that  submission,  the  paper  was  reviewed  by  two 
reviewers  and  I  am  in  receipt  of  their  reviews.  The  first  reviewer  admits  that  he  is  long- 
winded  and,  in  addition,  he  is  patronizing.  For  the  long-windedness,  he  is  forgiven.  I  am 
not  going  to  review  either  his  review  or  that  of  the  second  reviewer.  It  suffices  to  say,  in 
response,  that  the  issue  has  not  been  centered  about  whether  an  integral  can  or  cannot  be 
performed.  The  issue  is  more  subtle  than  that:  Has  the  integral  been  obtained  legitimately 
and  does  that  result  make  sense?  The  revised  manuscript,  enclosed  herein,  tries  to  make 
this  very  point  clearer  than  the  original  manuscript  did.  This  was  particularly  requested  by 
the  second  reviewer.  I  also  omitted,  for  now,  the  energetics  of  this  complex  structure. 
These  will  be  submitted,  in  due  course  under  separate  cover.  Thus,  it  may  be  possible  that 
the  reviewers  have  positively  contributed.  May  I,  nonetheless,  apologize  to  the  authors  of 
Reference  1  (now  Reference  (5)).  It  was  not  my  intention  to  torpedo  their  paper  or  imply 
that  they  are  wrong.  I  said  it  then,  and  I  repeat  it  now:  Once  the  auxiliary  condition  of 
modal  overlap  is  addressed  and  satisfied,  References  1-6  are  largely  validated.  I  do 
believe,  however,  that  this  condition  is  implicitly  made  in  all  these  references  and, 
therefore,  setting  the  loss  factors  of  the  satellite  oscillators  a  priori  to  zero,  is  a  violation  of 
what  has  already  been  implied. 


IV 


I  have  no  objection  to  this  enclosed  version  of  the  paper  being  published  as  a  letter, 
although,  I  would  prefer  a  paper,  subject  to  acceptance  in  the  first  place.  The  oversight  of 
the  auxiliary  modal  overlap  condition  has  been  circulating  about  for  more  than  five  years 
and,  in  spite  of  oral  pleas,  it  has  gone  unnoticed.  The  convulsion  over  “where  the  energy 
went”  is  still  an  issue  for  some.  It  is  time  to  lay  it  to  rest.  I  hope  this  paper  will  do  just 
that. 

Thank  you  for  your  interest  and  I  await  your  decision  in  this  matter. 

Sincerely  yours, 


G.  Maidanik 
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Abstract 


It  has  been  claimed  that  induced  damping  can  be  derived  from  the  coupling  of  a 
continuous  distribution  of  undamped  oscillators.  This  claim  stems  from  the  fact  that  the 
contribution  to  the  damping  of  a  master  oscillator  by  a  coupled  set  of  continuously  distributed 
satellite  oscillators,  is  independent  of  the  loss  factors  of  the  oscillators  in  this  set.  (The 
distribution  is  with  respect  to  the  frequency  of  resonance  of  the  satellite  oscillators  in  the  set). 
In  the  determination  of  the  induced  damping,  the  transition  from  a  discrete-to-a  continuous 
distribution,  however,  cannot  be  achieved  without  the  imposition  of  modal  overlap  on  the 
distribution  and  on  the  inherent  damping  of  the  satellite  oscillators.  It  is  this  imposition  that 
ensures  that  the  contribution  to  the  damping  by  these  satellite  oscillators  is  intuitively  real.  The 
imposition  forbids  equating  the  loss  factors  of  the  satellite  oscillators  to  zero  just  because  their 
contribution  to  the  surrogate  damping  —  the  induced  damping  --  of  the  host  master  oscillator  to 
which  they  are  coupled,  is  independent  of  these  loss  factors. 
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Introduction 


A  central  theme  to  a  few  publications  is  that  a  set  of  satellite  oscillators  attached  to  a 
master  oscillator  contributes  to  the  damping  of  the  master  oscillator  even  if  all  the  loss  factors  of 
the  satellite  oscillators  in  the  set,  are  equal  to  zero  [1-6]  [cf.  Fig.  1.]  The  satellite  oscillators 
are  numerable  and  their  resonance  frequencies  are  distributed  on  both  sides  of  the  resonance 
frequency  of  the  master  oscillator.  The  conclusion  that  a  set  of  lossless  satellite  oscillators  can, 
nonetheless,  contribute  to  the  damping  of  a  master  oscillator,  has  been  supported  by  analyses 
that  a  priori  assume  the  satellite  oscillators  to  posses  zero  loss  factors  [4, 5].  In  contrast  to  the 
analysis  presented  here,  some  of  the  analyses;  e.g.,  the  analysis  in  Reference  5,  are  conducted 
in  the  time  domain  [3].  With  respect  to  the  material  presented  in  this  paper,  whether  the 
analysis  is  conducted  in  the  time  domain  or  in  the  frequency  domain,  is  considered  to  be  merely 
a  matter  choice,  notwithstanding  that  in  the  case  of  a  closely  packed  distribution  of  resonance 
frequencies  of  the  satellite  oscillators,  the  relaxation  times  involved  may  be  quite  long.  In  this 
case,  to  relate  the  results  in  the  time  domain  to  those  in  the  frequency  domain  may  require  long 
times  to  capture  the  entire  data  [3]. 

In  this  present  paper,  a  fundamental  initial  oversight  that  beset  all  these  publications  is 
brought  to  the  attention  of  the  reader.  Admittedly,  once  that  oversight  is  addressed,  the 
remaining  analyses  and  arguments  in  all  these  publications  are  largely  validated;  e.g.,  validated 
is  the  argument  that  specifically  defined  measure  of  the  damping  that  is  provided  to  the  host 
master  oscillator  by  a  distribution  of  nearly  lossless  oscillators,  is  independent  of  the  loss 
factors  that  account  for  the  dampings  in  these  oscillators.  This  independence  does  not, 
however,  imply  that  individual  loss  factors  that  are  equal  to  zero  are  admissible;  indeed, 
individual  loss  factors  that  are  not  equal  to  zero  are  an  essential  condition  on  the  validity  of  the 
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analysis  in  the  first  place  [2,7].  Moreover,  the  analysis  of  nearly  continuous  distribution  of 
nearly  undamped  satellite  oscillators  brings  insights  into  the  manner  by  which  the  satellite 
oscillators  contribute  damping  to  the  mechanical  system  comprising  a  hosting  master  oscillator 
and  these  satellite  oscillators.  In  part,  these  insights  are  obscured  in  an  analysis  that  is 
subjected  a  priori  to  limiting  asymptotic  conditions  of  negligibility  and  continuity; 
notwithstanding  that  such  impositions  may  require  careful  consideration  if  apparent  singular 
behaviors  of  the  mechanical  system  are  to  be  avoided.  However,  detailed  considerations  of  this 
kind  are  to  be  presented  under  separate  cover  [8].  Here,  as  stated,  the  focus  is  on  the  linear 
equations  of  motion  that  describe  the  response  of  a  complex  structure  comprising  a  master 
oscillator  and  a  set  of  coupled  satellite  oscillators.  The  induced  loss  factor  in  the  damping  of 
the  master  oscillator  is  of  particular  interest.  It  is  argued  that  the  determination  of  the  induced 
loss  factor  is  readily  achieved  by  replacing  integration  for  the  summation  over  satellite 
oscillators.  However,  such  a  replacement  demands  the  obedience  of  an  auxiliary  relationship. 
This  relationship  imposes  that  a  modal  overlap  condition  be  satisfied.  A  satisfied  condition  of 
modal  overlap  states  that  if  the  loss  factors  of  the  satellite  oscillators  are  made  to  approach 
infinitesimal  values,  the  number  of  these  satellite  oscillators  must  approach  infinity  at  least  at 
the  same  rate.  This  condition,  then,  precludes  setting  the  loss  factors  of  the  satellite  oscillators 
equal  to  zero  a  priori  and  in  that  sense,  the  imposition  of  the  descriptive  adjective  “nearly”  is 
meaningful.  Conversely,  it  is  asserted  that  adhocly  replacing  a  summation  by  an  integration 
renders  the  condition  of  modal  overlap  implicitly  and  a  priori  satisfied.  The  significance  of 
these  statements  is  that  if  the  condition  of  modal  overlap  is  satisfied,  collectively  the  satellite 
oscillators  possess  enough  damping  to  explain  the  physical  presence  of  the  induced  loss  factor. 
Then,  inventing  exotic  mechanisms  to  account  for  the  induced  damping  becomes  moot. 
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I. 


The  Formalism  of  a  Discrete  and  of  a  Continuous  Distribution  of  Resonance 


Frequencies. 


Consider  the  mechanical  system  comprising  a  master  oscillator,  with  mass  (M)  and 
stiffness  ( K ),  that  is  coupled  to  a  set  of  satellite  oscillators.  This  mechanical  system  is 
sketched  in  Fig.  1  [1-6].  The  (r)th  satellite  oscillator  is  defined  by  a  mass  (Wlr)  and  a 
stiffness  (kr).  The  damping  is  assumed  to  be  associated  with  the  stiffness  elements 


K  =  K0(l  +  iTJ0)  ; 

(K0  /  M)  =  COq  , 

(la) 

kr  =  kor(l+iTJr)  ; 

{k0r  /  wiy)  =  (Of-  , 

(lb) 

where  (J]0)  and  (J]r)  are  loss  factors.  Each  of  these  loss  factors  characterizes  an  individual 
oscillator  and  quantifies  the  dissipation  that  the  oscillator,  in  isolation,  can  handle  when 
externally  driven.  The  linear  equations  of  motion  of  the  master  oscillator  in  situ  and  of  a  typical 
satellite  oscillator  in  situ  are 


R 

[(iCOM)  +  (K  /  i0))]Vom  +  %(kr/  i0))[Vom  ~  Vr(CO)]  =  Pe(G>)  .  (2) 

( iCdmr )  Vrm  +  (kr  /  iCO)  [Vrm-  Vo(0))]  =  0  ,  (3) 

respectively,  where  VQ((0)  and  Vr(CO)  are  the  responses  of  the  mass  (M)  of  the  master 
oscillator  and  of  the  mass  (mr)  of  the  (r)th  satellite  oscillator,  respectively,  (R)  is  the 
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number  of  satellite  oscillators  that  are  coupled  to  the  master  oscillator  and  P£  (CO)  is  the  drive 
that  is  assumed  applied  externally  to  the  master  oscillator;  the  satellite  oscillators  are  not  driven 
externally  [7].  [cf.  Fig.  1.]  The  satellite  oscillators,  in  Eqs.  (2)  and  (3),  are  assumed 
uncoupled  to  each  other;  the  satellite  oscillators  are  coupled  to  the  master  oscillator  only.  The 
summation  in  Eq.  (2)  is  over  a  set  of  oscillators,  where,  again,  (R)  is  the  number  of  satellite 
oscillators  in  the  set.  From  Eqs.  (2)  and  (3)  one  obtains 

Zo(C0)Vo(C0 )  =  Pe(CO)  ;  Vr(CO)  =  B(y,Xr,7]r)V0(CO)  ,  (4) 

where 

z0(co )  =  (iC0M){\-(yr2[i-s(y)]-i(yr2[rio+ris(y)]}  , 

R 

[s(y)-iris(y)}  =  (y)2  X  A(y,mr,xr,T]r) 

l 

A(y,mr,xr,T]r)  =  mrB(y,xr,T}r )  , 

B(y,xr,rjr)  =  (xr)2[(i  +  7 fr)(xr)2  -  (y)2  -  irir(y )2] 

{[(xr)2-(y)2]2+[lr(xr)2frl  .  (7) 

y  —  (CO  /  C0q)  ;  YYiy  —  (trip  /  M)  ;  Xp  =  (C0p  /  c Oq )  ;  1  ^  t  ^  R  .  (8) 

The  quantity  (Xp)  defines  the  resonance  frequency  distribution  of  the  satellite  oscillators.  It  is 
emphasized  that  in  the  present  paper  this  distribution,  and  not  the  spatial  distribution  of  the 


(5) 


(6) 
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satellite  oscillators,  is  in  reference.  The  satellite  oscillators  are  attached  on  a  single  velocity 
platform;  that  platform  constitutes  the  mass  of  the  master  oscillator.  The  velocity  in  the  plane  of 
the  platform  is  uniform  so  that  each  attached  sprung  mass  (satellite  oscillator)  perceives  the 
same  velocity  at  the  point  of  attachment.  Whether  the  spatial  distribution  of  the  sprung  masses 
on  the  platform  is  made  to  coincide  with  the  resonance  frequency  distribution  is  moot  in  this 
paper,  [cf.  Fig.  1.]  In  some  publications  variations  on  this  theme  are  implied  but  are  not 
always  strictly  defined.  Thus,  the  problem  of  a  priori  assuming  a  spatially  discrete  distribution 
of  satellite  oscillators  (sprung  masses)  versus  a  spatially  continuous  distribution  was  briefly 
discussed  in  Reference  8.  In  this  reference  the  satellite  oscillators  all  possess  the  same 
resonance  frequency.  In  Reference  6,  on  the  other  hand,  a  single  couple  structure  with 
multiple  resonance  frequencies  is  used. 

It  may  be  useful  to  exemplify  much  of  the  analytical  processes.  For  this  purpose,  two 
distinct  resonance  frequency  distributions  for  the  satellite  oscillators;  i.e.,  (Xf)  as  a  function  of 
(r),  are  depicted  in  Fig.  2.  In  Figs.  2a  and  2b  the  resonance  frequency  distributions  are  [5,  7] 

0 C0r/0)o )  =  Xr  =  [l+{l-2r}yCR)r1/2  \  \<r<R  ,  (9a) 


and 


0 C0r/C0o )  =  Xr  =  r(l-r)  1  ;  1  <r<R 


(9b) 


respectively,  where 


r  =  r(R+\Tx  ;  R  =  R(R+1T1  ;  y(R)  =  (7/2R)<(l/2)  .  (10a) 
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In  keeping  with  References  7  and  5,  a  normalized  mass  distributions  for  the  satellite  oscillators; 
i.e.,  (mr)  as  a  function  of  (r),  are  also  introduced  and  displayed  in  Fig.  2.  In  Figs.  2c  and 
2d  the  normalized  mass  distributions  are 

_  R 

mr  =  (MS/M)(R)-1  ;  (Ms/M)  =  £  mr  ,  (9c) 

1 

Mr  =  (Ms/M)[2/m  +  Vm-r)2+(r)2]-1  ; 

R  _ 

(Ms/M)  =  X  mr  ,  (9d) 

1 

respectively,  where  (Mg  /  M)  is  the  mass  ratio  of  the  total  mass  of  the  satellite  oscillators  to 
that  of  the  mass  of  the  master  oscillator.  In  this  paper  this  mass  ratio  is  set  at  one-tenth; 

(Mg/M)  =  (10  r1. 

The  impedance  Z0(CO),  stated  in  Eq.  (5),  comprises  the  self-impedance  ( i(OM ) 

[1 -  O0_2(l  +  °f  the  master  oscillator  and  the  sum  over  the  impedances  contributed  by 
the  satellite  oscillators  to  which  the  master  oscillator  is  coupled  [7].  Typically,  the  impedance 
C i(OM )  A(y,mr,xr,7jr)  is  that  contributed  by  the  (r)th  oscillator;  it  is  a  function  of  the 
normalized  frequency  y  =  (CO  /  C0o)  and  is  a  functional  of  the  normalized  mass  distribution 
ftlr  =  (rnr  /  M),  the  normalized  resonance  frequency  Xr  =  (COr  /  C0o)  and  the  loss  factor 
(T\r)  of  the  (v)th  satellite  oscillator.  As  Eqs.  (la)  and  (8)  state,  the  normalizing  frequency 
(C0o)  is  the  resonance  frequency  of  the  master  oscillator  in  isolation.  Under  the  same  cover, 
the  response  Vr(CO),  of  the  mass  (Wlr)  of  the  (r)th  satellite  oscillator,  is  related  to  the 
response  V0((D) ,  of  the  mass  (M)  of  the  master  oscillator,  by  a  factor  that  is  a  function  of  the 
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normalized  frequency  y  =  (CO  /  C0o)  and  is  a  functional  of  two  of  the  three  parameters; 
namely,  of  (Xr)  and  (Tjr).  The  relationship  between  Vr(CO)  and  V0(CQ)  is  stated  in  Eq.  (4). 

Clearly  (r)  is  a  discrete  variable.  Can  an  artificially  assigned  (r)  set  to  be  a  continuous 
variable?  This  proposed  continuity  of  (r)  is  exemplified  by  the  solid-line  curves  in  Fig.  2.  It 
is  convenient,  in  this  connection,  to  normalize  (scale)  the  continuous  variable  ( r )  by  the 
number  of  satellite  oscillators  (R)  plus  unity  (1),  as  suggested  by  Eqs.  (9)  and  (10).  The 
normalized  (r)  is  designated  (r )  and  the  dependent  parameters  are  then  expressed  in  the  forms 

mr=>m(F)  ;  xr=*x(r)  ;  T]r=$r\(r)  .  (11) 

The  scaled  forms  of  (Xr) ;  namely,  X(r),  in  Eqs.  (9a)  and  (9b),  are  depicted  as  functions  of 
(r)  in  Figs.  3a  and  3b,  respectively.  In  this  vein,  Eqs.  (4)-(7)  may  be  recast,  without  further 
ado,  in  the  forms 

z0((D)V0((D)  =  Pe((0)  ;  V(r,co)  =  B{y,  x(r),  mmco)  ,  (12) 


where 


z0((0)  =  (iC0M){i-(yy2[i-s(y)]-i(yr2[rio+ris(y)]}  (13) 

C(R+£)  __ 

[s(y)-iT]s(y)]  =  O02J(£)  dr  A{y,  m(r),  x(r),  7?(r)}  ; 

A{y ,  m(r),  x(r),  77(F)}  =  M(r )  B{y,  x(r),  7 7(F)}  ,  (14a) 
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C(R+£) 

[s(y)-iris(y)]  =  OO2.)^  dr  A{y,  H(r),  x(r),  7](r)}  ; 

A{y,  fl(r),  x(r),  77(F)}  =  }!{?)  B{y,  x(r),  77(F)}  ,  (14b) 

B{y,  x(r),  77(F)}  =  {x(F)}2  [(1+ {77(F)}2)  {x(r)}2  -(y)2  -iri(r)(y)2] 

{[{x(f)}2-(y)2f  +  [7 7(F){x(F)}2]2}1  ,  (15) 

mr  =*■  fl(r)  (i?  +  l)_1  ,  S  =  £(R  + 1)-1  ^  =  R(R  +  \)-'  .  (10b) 

The  parameter  (£)  designates  the  end-effects  that  are  introduced  in  the  transition  from  the 
discrete  to  the  continuous  domains;  e.g.,  £  —  £(R  + 1)-1 ,  with  £  ~  1/2,  say. 

[cf.  Fig.  2.] 

Afresh,  consider  the  mechanical  system  comprising  of  a  master  oscillator,  with  mass 
(M)  and  stiffness  ( K ) ,  that  is  coupled  to  a  set  of  satellite  oscillators,  as  defined  in  Eq.  (1)  and 
sketched  in  Fig.  1.  However,  now  the  resonance  frequencies  { &)(£)}  of  the  satellite 
oscillators  are  assumed  to  be  a  priori  continuously,  but  not  necessarily  uniformly,  distributed  in 
the  interval  0  <  £  <  1  [5].  Example  of  two  typical  resonance  frequency  distributions;  i.e., 

[JC (f)]  as  functions  of  (£),  are  depicted  in  Fig.  3.  In  Figs.  3a  and  3b  the  normalized 
resonance  frequency  distributions  are 

[co(0/(o0\  =  x(0  =  [i+(i-2 OY(R)Tm)  ;  (R)<(  1/2)  ;  o<c<i  ,  (16a) 
[(0(0/ co0]  =  x(0  =  (Od - 0“ 1  ;  o<c<i  ,  (16b) 
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respectively,  and  y(R)  is  defined  in  Eq.  (10)  [5, 7].  [cf.  Eqs.  (9a)  and  (9b),  respectively,  and 
Eq.  (10)].  The  corresponding  normalized  mass  distributions  /!(£)  of  the  satellite  oscillators 
are  [5  7,] 


H(0  =  (Ms  /  ;  0<f<l  ,  (16c) 

IKO  =  (Mslimin)  [(l-f)2+(f)V  ;  0<f<l  ,  (16d) 


[cf.  Eqs.  (9c)  and  (9d),  respectively,  and  Eq.  (10).]  The  linear  equations  of  motion  of  the 
master  oscillator  insitu  and  of  a  typical  satellite  oscillator  insitu  are 


z0(co )  vo(0))  =  pe(0) )  ;  v(C,co)  =  B{y,  x(0 ,  77(C)}  v0(®)  .  d?) 


where 


z0m  =  (io)M){i-(yy1li-sc(y)]-KyT2[rio+ris(y)])  ,  («) 

[Jc(y)-/J)f(y)]=0>)2  Jo  d;  A{y,  /«f),  x(f),  r/(0)  ; 

A{(y,  /1(f),  x(0,  >7(0)  =  W0«{y.  *(0,  ij(f))  ,  (19) 
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my,  x(a  mo)  =  wo)2  [(i+  n2( oiwof-w2  -mxyfi 


awo)2  -  w2]2+  tmowo)2)2)'1 .  (20) 

y  =  (colco0 )  ;  X(0  =  [OKO/Ool  ;  0<f<l  .  (21) 

A  question  arises:  Under  what  conditions  can  the  linear  Eqs.  (4)-(7)  be  asymptotically 
circumvented  into  the  linear  Eqs.  (12)-(15)  and  Eqs.  (12)-(15)  into  the  linear  Eqs.  (17)-(20), 
respectively,  and  are  these  conditions  significant?  This  question  has  been  partially  answered  in 
previous  publications;  in  this  paper,  a  more  definitive  answer  is  sought  [2,7].  However,  prior 
to  this  pursuit,  one  may  inquire:  Suppose  these  asymptotic  circumventions  were  found 
permissible  as  these  are  stated,  can  one  perform,  a  priori,  the  integration  in  the  continuous 
domain  that,  in  fact,  replaces  the  summation  in  the  discrete  domain? 
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II.  Evaluation  of  the  Induced  Loss  Factor  [7]s  (y)]  in  Eqs.  (13)-(15) 

The  following  quantities  are  defined  first 

z(F)  =  [x(r)/y];  [d{x(f)}ldrT'  =  f(f )  .  (22) 

Then,  from  Eqs.  (13)  -  (15),  one  obtains 

fZ(i?+£) 

=  W3Jz(£)  rfj(r)[/(F)/i(F)][{z(F))2l7(r)]  • 

{[{z(F))2-1]2+[(z(F)}27KF)]2)-‘  .  (23) 

Were  the  factor  [/(r)jti(r)]  to  be  a  well  behaved  function  of  (r)  in  the  range  £  <  r  <  R  +  £ 
and  were  7](r)  to  be  small  compare  with  unity,  as  a  function  of(r)  in  that  same  range,  the 
integral  in  Eq.  (23)  immediately  yields  the  result 

risiy)  =  (7T/2 )(y)3[f(r0)H(r0)\  ;  z(r0)  =  1  ;  (£)<r0<(R  +  €)  ,  (24) 

which,  clearly,  is  independent  of  TJ(r)  provided  these  loss  factors  are  set  small  compared  with 
unity  in  the  range  (£)<r  <R  +  £  [2].  [The  well  behaved  function  [/  r)]  of  (r)  and 

the  smallness  of  the  loss  factor  7 7(r)  as  compared  with  unity,  both  in  the  range 
(£)  <  r  <  (R  +  £)  are  specified  as  such  in  order  to  ensure  that  the  contribution  to  the  integral 
by  the  resonance  peaks  of  the  integral  at  Z(r)-\  are  dominant.  Were  contributions  by  the 
peaks,  in  certain  ranges  of  (y),  not  to  be  dominant,  the  integral,  as  just  performed  will 
naturally  underestimate  the  value  of  the  integral  in  these  ranges  of  (y) .  In  the  present  paper, 
situations  that  may  fall  into  this  category  are  neither  evaluated  nor  specifically  considered. 
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Indeed,  this  is  the  reason  that  only  the  induced  damping  is  evaluated  and  considered  in  this 
paper  [8].]  Employing  the  normalized  resonance  frequency  distributions  stated  in  Eqs.  (9a) 
and  (9b),  one  finds  from  Eqs.  (22)  and  (24)  that 

f(f0)  =  [r(*)W3r‘  ;  [i+(r/2)r<1,2><:y<[i-(r/2)r<1,2>  ,  asu 

f(r0 )  =  (1  +  yr2  ;  (2R  + 1)-1  <y<(2R  +  l)  ,  (25b) 

respectively,  where  (y)  is  defined  in  Eq.  (10)  and  is  chosen,  herein,  to  be  equal  to  (0.6). 
Similarly,  employing  the  normalized  mass  distributions  stated  in  Eqs.  (9c)  and  (9d),  one  finds 
from  Eqs.  (22)  and  (23)  that 


IMT0)=*(MSIM){R)-X  ;  (£)  <  rQ  <(*  +  £),  (25c) 

LL(r0)=>(Ms/M)(2/Jl)(\  +  y)2  [l  +  OO2]"1  ;(£)  <  r0  <  (R+£)  .  (25d) 

From  Eqs.  (24)  and  (25)  one  readily  derives,  for  the  two  examples  herein  considered,  the 
results 

risiy)  =  (Ms  I  M)[K  /  {2nm  ;  [l  +  (y/2)r(1/2)<  y<  [l-(y/2)r(1/2)  ,  (26a) 

ris(y)  =  (M5/M)(y)3[l  +  (y)2r1  ;  [(2R+1)]-'  <  y  <  [(2i?+i)]  .  (26b) 

It  is  recognized  that  Tjs(y)  is  largely  meaningful  only  at  and  in  the  vicinity  of  (y)  equal  to 
unity;  y  ~  1,  so  that  the  restricted  range  of  (y)  in  Eq.  (26)  is  not  physically  significant.  For 
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the  two  examples,  stated  in  Eqs.  (26a)  and  (26b),  T]s(y)  is  displayed  in  Figs.  4a  and  4b, 
respectively.  The  dashed-line  curves  in  these  figures  are  merely  free  extrapolations  beyond  the 
permissible  ranges. 

Since  it  emerges  that  Tfs(y)  is  independent  of  the  loss  factors  Tj(r),  provided 
7]{r) « 1  in  the  range  (£)  <r  <(R  +  £),  these  loss  factors  of  the  satellite  oscillators  may, 
indeed,  be  set  equal  to  zero  and  nobody  would  be  the  wiser.  Is  that  true? 

Finally,  it  is  remarked  that  the  difference  between  the  induced  loss  factor  7js(y),  stated 
in  Eq.  (14),  and  Tl§(y)  stated  in  Eq.  (19),  is  the  number  (R)  of  satellite  oscillators.  Clearly, 
Eq.  (19)  can  be  derived  from  Eq.  (14)  by  allowing  ( R )  to  increase  enough  to  render 
(/?+£)->!  and  (£)^0. 
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III.  Modal  Overlap  Condition  for  Relating  the  Discrete  and  the  Continuous  Domains 

To  establish  the  relationship  between  the  discrete  and  the  continuous  domains  one  may 
define  two  distinct  frequency  bands.  The  first,  designated  (8(X)r)  defines  the  fair  frequency 
territory  occupied  by  each  satellite  oscillator  with  respect  to  its  adjacent  neighbors;  this 
frequency  bandwidth  is  centered  on  the  resonance  frequency  of  that  satellite  oscillator.  The 
expression  for  this  bandwidth  is 

SC0r  =  {(l/2)(fi)r+1-G)r_1)  [U(R-l  +  £-r)-U(2-£-r)] 

+  (fi)2-ft)i)t/(l  +  e-r)  +  (COR-COR-iMr+e-R^Uir-SMR+e-r)  ,  (27a) 

where  U(OC)  is  the  usual  unit  step  function  and  (£)  is  less  than  unity,  equal  to  (1  /  2),  again, 
say.  The  inverse  of  the  frequency  bandwidth  ( 8(Or )  is  commonly  referred  to  as  the  modal 
density  ( Ylr )  of  the  satellite  oscillators;  i.e., 

nr  =  (8(Or)~l  ;  e<r<R+£  ,  (28) 

where  the  inversion  of  ( 8(Or )  is  carefully  executed  with  respect  to  the  unit  step  functions  [10, 
11].  The  second  frequency  band,  designated  (A COr),  defines  the  inherent  frequency 
bandwidth  of  the  (r)th  satellite  oscillator;  this  bandwidth  accounts  for  the  inherent  damping  of 
this  satellite  oscillator.  The  measure  of  the  damping,  in  turn,  defines  the  loss  factor  (T}r)-  The 
bandwidth  (A CQr)  is  then  stated  in  terms  of  the  loss  factor  (Tjr)  and  the  resonance  frequency 
(C0r)  in  the  form 
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ACOr=(CQrrir)U(r-£)U(R+£-r)  . 


(27b) 


The  ratio  (br)  of  these  two  frequency  bandwidths  defines  a  modal  overlap  parameter 

{br  =  (A COr/SCOr)  =  (I nrCDrTlr)}U(r-£)U(R  +  £-r )  ,  (29) 

where  use  is  made  of  Eqs.  (27)  and  (28). 

For  a  modal  overlap  parameter  (br)  that  is  less  than  unity;  br<  1,  adjacent  satellite 
oscillators  reside  outside  each  other’s  bandwidths.  Consequently,  the  influence  of  the  satellite 
oscillators  on  the  response  of  the  master  oscillator,  as  a  function  of  (y);  y  =  (0)  /  COq),  can  be 
identified  individually;  each  contribution  associated  with  a  satellite  oscillator  stand  out 
prominently  from  the  others.  An  example  of  such  an  influence  in  terms  of  evaluating  the 
induced  loss  factor  Tl$(y),  as  a  function  of  (y)  with  (br)  less  than  unity,  is  depicted  in 
Fig.  5  [5, 7].  On  the  other  hand,  for  a  modal  overlap  parameter  (br)  that  exceeds  the  value  of 
unity;  br>  1,  adjacent  satellite  oscillators  reside  within  each  other’s  bandwidths.  Therefore, 
their  influence  on  the  response  of  the  master  oscillator  is  largely  continuous  as  a  function  of 
(y)-,  y-(CO/  C0o).  The  more  the  modal  overlap  parameter  (br)  exceeds  the  value  of  unity, 
the  more  the  continuity  [7].  An  example  of  such  an  influence,  in  terms  of  evaluating  the 
induced  loss  factor  Tjs(y),  as  a  function  of  (j)  with  (br)  exceeding  unity,  is  depicted  in 
Fig.  6  [5,  7].  Comparing  Fig.  4  with  Fig.  5  and  with  Fig.  6,  respectively,  serves  to  expose 
the  issue  that  is  central  to  this  paper.  Figures  4a  and  4b,  respectively,  overlap  prime  portions 
of  Figs.  6a  and  6b  only.  These  prime  portions  are  defined  within  the  range  set  on  (y)  in 
Eqs.  (26a)  and  (26b),  respectively.  No  such  overlap  exists  between  Figs.  4a  and  4b  and 
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Figs.  5a  and  5b,  respectively.  Nonetheless,  if  the  undulations  in  Fig.  5  are  appropriately 
averaged,  the  overlap  between  Fig.  4  and  6  is  extended  to  include  the  mean  values  of  Fig.  5 
[7,  12].  [cf.  Appendix  A.]  When  (br)  becomes  small  compared  with  unity  the  mean  values 
alone  do  not  reflect  the  presence  of  undulations  and,  therefore.  Fig.  4  cannot  serve  to 
substitute  for  Fig.  5.  The  information  that  lies  in  Fig.  5  can  hardly  be  derived  from  Fig.  4; 
Fig.  4  is  akin  to  Fig.  6  but  not  to  Fig.  5.  Since  Fig.  4  is  derived  by  replacing  the  summation 
by  integration,  it  is  concluded  that  this  replacement  is  commensurate  with  the  imposition  of 
modal  overlap  parameters  that  exceed  unity;  by>  1.  [cf.  Appendix  A.]  An  imposition  of  this 
kind  ensures  that  the  evaluations  of  mean- value  data  do  not  conceal  information  that  may  render 
the  data,  at  best,  ambiguous  and,  at  worst,  misleading  [8, 13]. 

When  the  modal  overlap  parameter  (br)  exceeds  unity,  the  modal  overlap  condition  is 
said  to  be  satisfied.  When  the  modal  overlap  condition  is  satisfied,  the  linear  dynamic 
description  [S()0  -  iTJgiy)] ,  comprising  the  sum  of  the  individual  linear  dynamic  descriptions 
of  all  the  satellite  oscillators;  e.g.,  typically  the  term  [(y)2  A(y,mr,Xf,T]f)]  contributed  to  the 
sum  by  the  ( r)th  satellite  oscillator,  may  be,  well  nigh,  evaluated  by  integration 

R  _ 

[s(y)-iris(y)]  =  OO2  X  A(y,mr,xr,rir) 

l 

(•(*+£)  _ 

=>  (y)2J(£)  drA{y,m(r),x(r),ri(r)}  ,  (30a) 

where  (/*),  in  the  second  of  Eq.  (30a),  is  a  continuous  and  dimensionless  variable  [2]. 

[cf.  Eqs.  (6)  and  (14a).]  It  is  emphasized  that  the  transition,  described  in  Eq.  (30a),  is 
validated  if,  and  only  if,  it  is  implicitly  understood  that  (br),  as  stated  in  Eq.  (29),  exceeds 
unity;  namely 
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{ br  =  {nrCOrT]r )  >  1}  U(r-£ )  U(R  +  £-r )  . 


(31a) 


The  assignment  of  continuity  to  (r),  in  Eq.  (30a),  renders  the  normalized  resonance  frequency 
Xr  =»  X(r)  a  continuous  function  of  (r),  as  exemplified  in  Figs.  2a  and  2b  [14].  Similarly, 
mr  =>  m(r)  and  Tjr  =>  Tj(r).  [cf.  Figs.  2c  and  2d.]  In  this  rendering,  Eqs.  (27a)  and 
(27b)  may  be  approximated,  by  smoothing,  in  the  forms 

n(r)  =  [ftOWr1  =  [< dCOir )  /  drY1  U(r  -£)  U(R+£-r)  ,  (32a) 

AO)(r)  =  (D(r)7](r)  ;  (£)<  r<(R  +  £)  ,  (32b) 

respectively,  where  one  recognizes  that  [ dO)(r )  /  d(r)]  is,  by  definition,  a  positive  quantity 

and  Eq.  (28)  is  referenced.  From  Eq.  (22)  one  finds  that 

[dCQ(r)  /  d(r)]  =  [(/?  +  l)/(r)  /  (0oTl  ■  Under  the  guidance  of  Eq.  (32a)  and  the 

normalization  (scaling)  proposed  in  Eq.  (10),  Eq.  (9)  yields  the  two  corresponding  modal 

densities 


n(r)  =  n{r)  =  [(/?  +  l)/©07(i?)]{x(r)}"3  U(r-£ )  U(R  +  £-r )  ,  (33a) 

n{r )  =  n(r )  =  [(tf+l)/C%](l-r)2  U(r-£ )  U(R  +  £-r )  ,  (33b) 

respectively,  where,  again,  £  =  £(R  + 1)-1 ,  R  =  R(R  + 1)-1 ,  r  =  r(R  + 1)-1  and 
y(R)  =  (/  /  2 R) .  From  Eq.  (31a),  the  modal  overlap  condition  is  satisfied,  in  the 
continuous  domain,  provided 
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{b(r)=b(r)=n(r)CO(r)V(r)=n(r)CQ(r)ri(r)>l}  U(r-£ )  U(R  +  £-r )  .  (31b) 

In  terms  of  Eqs.  (33a)  and  (33b),  the  specific  modal  overlap  conditions,  per  Eq.  (31b),  are 

{b(r)=b(r)=[X(r)T2[(R  +  l)ri(F)/y(R)]>l  }  U(r-£ )  U(R  +  £-r )  ,  (34a) 

{b(r)=b(r)=[r(l-rMR+l)T](r)]>l}  U(r-£)  U(R+£-r)  ,  (34b) 

respectively.  It  is  emphasized  that  the  transition,  from  the  discrete  to  the  continuous  domains, 
state  in  Eq.  (30a),  is  predicated  on  the  validity  of  Eq.  (31),  in  general,  and  on  Eq.  (34) 
specifically  for  the  two  examples  here  pursued.  To  facilitate  the  assessment  of  the  transition, 
from  the  discrete  to  the  continuous  domains,  it  is  convenient  to  scale  Eq.  (30a).  The 
continuous  (r)  is,  in  that  scaling,  normalized  by  the  number  (R  + 1),  as  already  performed  for 
example  in  Eq.  (33).  The  so  normalized  continuous  (r)  is  designated  (£)  and  hence,  in  that 
vein,  Eqs.  (31)-(34)  are  to  be  recast  with  r  — >  (£);  e.g.,  in  this  normalization  Eq.  (32)  is 
written  in  the  form 

{n(OdO)(0  =  (R  +  l)dC)  U(£-£ )  U{R+£-Q  ,  (35a) 

ACO(C)  =  coiCMO  ,(£)<£<  ( R+£ )  ,  (35b) 

respectively,  where 


19 


(36a) 


£  =  r(/?  +  l)_1  ;  £  =  £  (R  +  i)'1  ;  R  =  R  (R+ 1)"1 

The  scaling  is  illustrated  in  Fig.  3  for  the  examples  depicted  in  Figs.  2a  and  2b.  [cf.  Eqs.  (9) 
and  (16).]  Employing  this  scaling  in  Eq.  (30a),  one  obtains 


f(K+£) 

[S(y)  - irjs(y )]  =  O02J(g)  dr  A{y,  m(r),  x(r),  7](r)} 

f(/?+£) 

=>  (yf  J(g)  dy %  A{y,  X(0,  m)}  .  (30b) 

where  (£)  and  (£)  are  defined  in  Eq.  (36a).  [cf.  Eqs.  (10)  and  (14).]  If  one  further 
assumes  that  the  number  (R)  of  the  satellite  oscillators  is  large  compared  with  unity,  Eq.  (30b) 
becomes 

[S09  -irjs(y)]  =>  (y)2J0d(  My*  140,  x(0,  7/(0) 

=  [ic(>>)  -ir}§(y)]  ,  (30c) 

where  [5C(j)  -  iTfs(y)]  is  defined  in  Eq.  (19)  and  it  follows  from  Eq.  (36a)  that 

£— >o  ;  (R  +  £)—»!  ;  R— >1  as  R  »>  1  .  (36b) 
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The  modal  overlap  condition  that  validates  Eqs.  (30b)  and  (30c)  is  derived,  from  Eq.  (31b),  in 
the  form 

{6(0  =  GKO  [d(0(C)/dcrl[(R+im)]>i }  utf-e)  u(R+e-0  ,  (3ic) 

where  use  is  made  of  Eq.  (35a).  For  the  examples  quoted  in  Eq.  (9),  Eq.  (31c)  assumes  the 
forms 

{b(0  =  wor2  K*+i)TO)/r(*)]>i}  U(;-S)  uOf+e-o  ,  (37a) 

(b(0  =  Kd-O)  [(«+ i)i7(0] >i)  V(i-£)  u(R+e-0  ,  (37b) 

where  in  Eq.  (37a) 


[x(or2  =  [i+a-2 orm  ;  ym  =  (y/2S)<a/2) 


(38) 


[cf.  Eq.  (16a).]  It  emerges,  therefore,  that  the  validity  of  Eqs.  (30b)  and  (30c)  demands  that 

Eq.  (31c),  m  general,  and  Eqs.  (37a),  and  (37b)  in  particular,  are  obeyed.  These  latter  three 
equations  may  be  cast  in  the  form 


{ K«+D  mo]> 


\[d\n{CQ(0}ldQ 

[7(*){*(0}2] 

Kxi-or1 


\  }  m-£)  U(R  +  £-0  , 


(3  Id) 
(37c) 
(37d) 
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The  quantity  [d  In (®(f) }  /  dCl  -  [d ln«0 )  /  df],  in  Eq.  (3  Id),  is  the  loeal  slope  in 
figures  of  which  Figs.  3a  and  3b  serve  as  two  specific  examples;  i.e.,  the  examples  given  in 
Eqs.  (37c)  and  (37d),  respectively  [14].  [It  is  noted  that  Eq.  (3 Id)  may  be  cast  in  the  alternate 

form 


7](r)  >  [d\n{CO(r)}/dr]  =  [d\n{x(r)}/dr]  .  (31e) 

In  Eq.  (3  le)  the  local  slope  is  that  of  curves  like  those  exemplified  in  Figs.  2a  and  2b  for  the 
continuous  (r)  [14].  Again,  it  is  emphasized  that  the  relationship  between  Figs.  2a  and  2b 
and  Figs.  3a  and  3b,  respectively,  is  merely  the  scale  factor  (*  +  D-  In  that  sense  Eq.  (3 Id)  is 
merely  the  normalized  version  of  Eq.  (31e)  and,  as  such,  the  former  is  to  be  preferred  when 
comparing  situations  of  varying  numbers  of  satellite  oscillators  [5].]  Equations  (37c)  and  (37d) 
are  depicted  graphically  in  Figs.  7a  and  7b,  respectively.  These  equations  and  figures 
exemplify  that  [77(01  cannot  be  independently  set  equal  to  zero  without  qualification.  The 
lower  limit  that  can  be  assigned  to  [77(01  is  dependent  on  the  number  (R)  of  satellite 
oscillators  plus  unity  (1) .  The  value  of  [<*  + 1)77(01  may  be  large  or  small  compared  with 
unity,  but  it  must  remain  finite  if  Eqs.  (30b)  and  (30c)  are  to  be  validated.  The  only  way  to 
allow  for  infinitesimal  [rj(0]’s;  i-e”  nearly  undamped  satellite  oscillators,  is  to  explicitly  or 
implicitly  render  (R  + 1)-1  proportionately  infinitesimal.  In  particular,  the  number  of  satellite 
oscillators  needs  to  approach  infinity  at  least  at  the  same  rate  that  [77(0]  approaches 
infinitesimal  values.  Thus,  once  the  modal  overlap  condition  is  satisfied,  the  satellite  oscillators 
collectively  possess  enough  damping  to  account  for  the  physical  presence  of  the  induced  loss 


factor. 
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Appendix  A:  A  Few  Additional  Remarks 


In  a  steady  state  analysis  of  the  dynamic  characteristics  of  the  response  of  structural 
complexes,  many  data  sets  are  presented  in  terms  of  appropriately  averaged  values.  Such 
presentations  are  ultimately  acceptable  as  long  as  these  mean  values  are  all  the  information  that 
is  either  available,  sought  or  both.  The  loss  factor  TJsiy)  presented  in  Fig.  6,  which  pertains 
to  a  modal  overlap  parameter  ( by )  that  is  equal  to  two;  by  -  2,  simultaneously  fits  the  two 
categories.  On  the  other  hand,  the  loss  factor  Tjs(y)  depicted  in  Fig.  5,  which  pertains  to  a 
modal  overlap  parameter  (by)  that  is  equal  to  one  tenth;  by  =  10~* ,  does  not  meet  either 
categories.  Figure  4,  however,  cast  the  data  in  Fig.  5  in  terms  of  mean  values,  thus  meeting 
the  availability  category  only;  the  data  sought  may  demand  some  form  of  accounting  for  the 
undulations.  The  undulations  are  excursions  in  the  values  of  TJs(y),  as  a  function  of  the 
normalized  frequency  ( y ) ,  between  resonance  peaks  and  anti-resonance  nadirs.  Since  the  size 
of  the  excursions  are  determined  by  the  modal  overlap  parameter  (by) ,  can  a  measure  of  the 
undulations  be  devised  as  a  function  of  this  parameter?  The  excursions,  it  is  noted,  are  the 
larger  the  smaller  are  the  modal  overlap  parameters  compared  with  unity.  One  is  reminded  in 
this  connection,  however,  that  to  a  first  order  of  approximation  results  obtained  from  averaging 
the  undulations;  e.g.,  the  appropriately  averaged  values  in  Fig.  5,  as  well  as  the  results 
obtained  by  performing  the  integration  in  Eq.  (23);  e.g.,  the  data  depicted  in  Fig.  4,  are  found 
to  be  independent  of  the  individual  loss  factors  of  the  satellite  oscillators  [1-7].  Therefore,  to 
that  approximation  the  mean  values  also  lack  dependence  on  the  modal  overlap  parameter  (br)\ 
e.g.,  although  there  is  no  demand  for  equal  values  of  this  parameter,  there  is  an  overlap 
between  Figs.  4  and  6.  On  the  other  hand,  although  the  values  of  the  modal  overlap 
parameters  may  be  identical  in  Figs.  4  and  5,  there  is  no  overlap  between  prime  regions  in 
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these  figures,  except  in  a  qualified  sense.  Indeed,  when  mean  values  only  are  to  be  reported,  to 
avoid  “a  tail  wagging  a  dog,”  it  is  to  be  assumed  that  the  modal  overlap  parameters  are  set  in 
excess  of  unity,  unless  they  are  otherwise  specifically  stated  [13].  Again,  when  values  are 
derived  by  appropriately  averaging  undulating  data,  the  value  of  the  modal  overlap  parameter 
(br)  must  accompany  the  specification  of  these  data;  e.g.,  when  Fig.  4  is  meant  to  represent 
the  mean  values  of  Fig.  5,  the  modal  overlap  parameter  ( by )  must  be  stated  as  one  tenth; 
br  - 10-1 .  Otherwise,  Fig.  4  is  meant  to  represent  Fig.  6  (for  which  by-  2)  and  not,  at  all, 
Fig.  5. 

In  this  connection,  it  transpires  that  the  first  order  approximation  just  discussed  is  the 
better  the  larger  is  the  number  (R)  of  the  satellite  oscillators  and  the  closer  is  the  modal  overlap 
parameter  (by)  to  unity.  Indeed,  for  low  values  of  (R)  and  for  large  values  of  (by),  erosions 
may  beset  the  exact  data  obtained  via  the  summation  in  Eq.  (6)  [7].  Such  simultaneous 
decreases  in  the  number  of  satellite  oscillators  and  increases  in  the  modal  overlap  parameters 
may  not  only  render  the  exact  data  dependent  on  these  quantities,  but  may  also  violate  the 
requirement  that  the  contributions  to  the  integral  be  dominated  by  the  resonances  in  the 
integrand,  [cf.  Eqs.  (23)  and  (24).]  Reference  7  suggests  that  to  ensure  a  first  order 
approximation  that  is  free  of  erosions  one  must  require  that  by  =  [br(R  + 1)-1  ] «  1 .  On  the 
other  hand,  Eqs.  (3  Id),  (37c)  and  (37d)  suggest  that  to  ensure  that  mean  values  are  the  prime 
data,  (by)  must  exceed  unity;  by  >  1. 
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Fig.  1.  Master  oscillator  coupled  to  a  set  of  resonance  frequency  distributed  satellite 
oscillators.  The  master  oscillator  is  defined  by  the  mass  (M)  and  the  stiffness 
(K)  and  the  (r)th  satellite  oscillator  by  mass  (1Tlr) and  stiffness  (kr).  Only 
the  master  oscillator  is  driven,  by  an  external  drive  Pe(CO),  generating  the 
response  Vq(CO)  in  the  master  oscillator  and  the  response  Vr(CO)  in  the  (r)th 
satellite  oscillator. 
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Fig.  2.  Resonance  frequency  distribution;  (Xj-),  and  normalized  mass  distribution, 

(mr),  as  functions  of  the  discrete  (r)  (distinct  dots),  and  X(r)  and  m(r),  as 
functions  of  the  continuous  (r)  (solid  line  curves).  Number  (R)  of  satellite 
oscillators  is  twenty  seven;  (R  =  27)  and  mass  ratio  (M$  /  M)  is  a  tenth  (10  ) 
a.  (Jtr)  as  specified  in  Eq.  (9a).  b.  (Xr)  as  specified  in  Eq.  (9b). 


m{r) 

d.) 

MO 

Fig.  2.  Resonance  frequency  distribution;  (Xr) ,  and  normalized  mass  distribution; 

(Mr),  as  functions  of  the  discrete  (r)  (distinct  dots),  and  X(r)  and  M(r),  as 
functions  of  the  continuous  (r)  (solid  line  curves).  Number  (R)  of  satellite 
oscillators  is  twenty  seven;  (R  =  27)  and  mass  ratio  (My  /  M)  is  a  tenth  (10-1) 
c.  (Wlr)  as  specified  in  Eq.  (9c).  d.  (m r)  as  specified  in  Eq.  (9d) 


normalized  continuous  variable  (r ) 
a.  Jt(£)  as  specified  in  Eq.  (16a). 
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Fig.  4.  Induced  loss  factor  Tjs(y ),  as  a  function  of  the  normalized  frequency  (J), 

resulting  from  replacing,  adhocly,  a  summation  by  integration  [/?  =  27;  (Ms  /  M)  =  10-1.] 

a.  For  a  normalized  resonance  frequency  {Xr)  and  a  normalized  mass 
distribution  ( mr )  as  specified  in  Eqs.  (9a)  and  (9c),  respectively. 

b.  For  a  normalized  resonance  frequency  (Xr)  and  a  normalized  mass 
distribution  ( Tnr )  as  specified  in  Eqs.  (9b)  and  (9d),  respectively. 
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Fig.  5.  Induced  loss  factor  Tfs(y) ,  as  a  function  of  the  normalized  frequency  ( J),  for  a 

modal  overlap  parameter  (br)  that  is  less  than  unity;  br  ~  (10)'1 .  [R = 27;  ( Ms  /  M) 

a.  For  a  normalized  resonance  frequency  (X/*)  and  a  normalized  mass 

distribution  (mr)  as  specified  in  Eqs.  (9a)  and  (9c),  respectively. 

b.  For  a  normalized  resonance  frequency  (Xr)  and  a  normalized  mass 
distribution  (Fnr)  as  specified  in  Eqs.  (9b)  and  (9d),  respectively. 


Fig.  6.  Induced  loss  factor  Tjs(y) ,  as  a  function  of  the  normalized  frequency  (y) ,  for  a 

modal  overlap  parameter  (by)  that  exceeds  unity;  by  ~  2.  [/?  =  27;  (My  /  M)  =  10  *.] 

a.  For  a  normalized  resonance  frequency  ( Xr )  and  a  normalized  mass 
distribution  ( mr )  as  specified  in  Eqs.  (9a)  and  (9c),  respectively. 

b.  For  a  normalized  resonance  frequency  ( Xr )  and  a  normalized  mass 
distribution  ( ntr )  as  specified  in  Eqs.  (9b)  and  (9d),  respectively. 


33 


Fig.  7.  Lower  limit  imposed  on  the  value  of  the  quantity  [(/?  + 1  )f/(£)]  >  35  a  function 
of  (£),  in  order  to  just  satisfy  the  modal  overlap  condition;  i.e.,  br  ~  l,  thereby, 
just  permitting  for  an  integration  to  replace  a  summation,  [cf.  Figs.  5  and  6.] 

(R  =  27) . 

a.  For  a  resonance  frequency  distribution  as  stated  in  Eq.  (16a). 

b.  For  a  resonance  frequency  distribution  as  stated  in  Eq.  (16b). 
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